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Abstract 

Symmetry breaking solutions are investigated in the N — > oo limit for the ground 
state of a system consisting of a Lorentz-scalar, N component "phantom" field and 
an O(N) singlet. The most general form of 0(N) x Z2 invariant quartic interaction 
is considered. The non-perturbatively renormalised solution demonstrates the pos- 
sibility for Z2 symmetry breaking induced by phantom fluctuations. It becomes also 
evident that the strength of the "internal" dynamics of the N-component field tunes 
away the ratio of the Higgs condensate and the Higgs mass from its perturbative 
(nearly tree-level) expression. 
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1 Introduction 



Scalar fields are proposed at present with renewed intensity as minimal, phenomeno- 
logically motivated complements to the Standard Model. The simplest mechanism for 
electroweak symmetry breaking is realised by introducing a renormalisable biquadratic 
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coupling of the Higgs field to a "hidden scalar" (the phantom field) which acquires nonzero 
vacuum expectation value [1,2]. Postulation of two extra scalars with appropriate dynam- 
ics might provide the minimal framework for the particle physics interpretation of the 
cosmologically well-established facts of missing luminous mass and of inflationary density 
fluctuations [3,4]. A scalar inflaton coupled to sterile right-handed neutrinos offers another 
cosmologically motivated minimal extension of the standard model, which accounts also 
for neutrino oscillation phenomena and the fermion- asymmetry of the universe [5]. 

A general background to this evolution is provided by low energy field theories constructed 
on string vacua showing maximal compatibility with the particle content of the Standard 
Model. Appearance of extra scalars is characteristic for all these constructions. A con- 
sequence of possible physical significance is that the natural cut-off scale to be applied 
for the top quark contribution to the Higgs mass would be raised considerably [6]. Con- 
sequences of the existence of a phantom sector on the precision electroweak tests of the 
Standard Model are being systematically explored [7]. The history of such investigations 
goes back to the 80's [8]. 

The goal of this letter is to illustrate the substantial influence of an extended scalar 
sector on the symmetry breaking vacuum features through the results of an exact, non- 
perturbatively renormalised large N analysis. We propose the following embedding of 
the Higgs sector of the Standard Model into a wider set of elementary scalar fields. The 
gauged self-interacting scalar SU (2) doublet is coupled also to a set of scalar fields which 
form an iV-component vector under a hidden internal O(N) symmetry. The coupling 
preserves the 577(2) symmetry of the standard gauge model. After fixing 3 components of 
the complex Higgs doublet to zero a Z2 invariance of the Higgs field is left. The phantom 
O(N) vector field has some extra O(N) invariant non-linear interaction on its own. With 
appropriately chosen TV-scaling of the quartic couplings the renormalised N = 00 solutions 
will be presented. We shall demonstrate that the non-perturbative renormalisation of the 
solution leads to a non-trivial Z(2) symmetry breaking Higgs condensate in a large range 
of the renormalised couplings of the model even if no negative renormalised squared mass 
is introduced. It will be shown that also essential deviations might appear in the relation 
of the Higgs mass to the strength of the Higgs condensate, relative to the leading order 
large N estimate derived in the O(N) symmetric linear sigma model. 

The present model is a generalisation of the 0(N + M) symmetric model studied by 
Chivukula et al. in 1991-1993 [9]. Another extension of the Higgs sector was investigated 
in [10], where consequences of the mirror world assumption were studied. Here the general 
model will be solved by taking the large TV limit of the renormalised Dyson-Schwinger 
equations for the exact n-point functions. The renormalisation process will follow an 
approach advocated recently in the context of 2PI approximate solution of scalar field 
theories [11], but we emphasise that it represents a rather non-trivial implementation of 
the strategy of nonperturbative renormalisation. 
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2 Leading order formal solution 



The Lagrangian to be solved at N = oo describes the most general renormalisable SU (2) x 
O(N) symmetric interaction of an iV-component (and SU(2) singlet) field ipi, i — 1, ..,N 
with a gauged SU(2) doublet (and O(N) singlet) $: 



Xl i^) 2 -^-r(^) 2 -^-M^- (1) 



2AN KTlJ 2AN K ' \2N 

The present parametrisation is a generalisation of that used in [9], where only a unique 
quartic coupling ~ (ip 2 + $^$) 2 was considered. 

We shall investigate the usual symmetry breaking pattern, where the lower real part 
of the complex SU{2) doublet gets a nonzero classical value y/~Nv, and the other three 
components are gauge transformed away (unitary gauge). The shift in the Higgs field then 
breaks the Zi symmetry related to the remaining reflection symmetry of the Higgs field. 
When gauge fields are set to zero the resulting shifted Lagrangian density is the following: 



Lfree[^, <A = \{d^f + - \m\$ - X -m\{cj 2 + 2ay/Nv + Nv 2 ), 

LUA, *] = -^ 2 f - ^(- 2 + *°VNv + Nv 2 ) 2 

A;i i/>?(o 3 + 2ay/Nv + Nv 2 ). (2) 
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The large N solutions will be constructed with help of the Dyson-Schwinger (DS) equa- 
tions which are generated by the master equation 



5F SS. 
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(3) 



and its functional derivatives [12]. The right hand side is evaluated by replacing in the ex- 
pression of the functional derivative of the classical action the field ipu by ipu + Gjjv8/8 1 Pv 
and the whole expression is applied to the unity in the field space. T is the effective action, 
ipu is a generic field variable, Gab is the two-point function of the generic fields <£a and 
(fB (deWitt's convention is used for the index A,B,...). Note, that one should include 
into the intermediate steps of the derivation all possible (A, B) pairs of indices for Gab, 
the corresponding external currents are set to zero only at the very end. 

Below we give the unrenormalised equations which determine the vacuum condensate 
v, the propagators G aa (x,y) =: G a (x,y), G Mn (x,y) =: 5 nm G^,(x,y) and the relevant 
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3-point function T^ n ^ ma (x,y, z) = 5 nm T^ a (x,y, z). The LO (N = oo) equation for the 
vacuum expectation value v looks rather simple: 

= 0, T^p = f G f (k). (4) 

Jk 

The LO tp propagator is obtained as 5 2 Y / 8ip8ip =: iG^ ) 1 {p) in the following form 

iG?{p) = v 2 - ml - - ±T,. (5) 

Since the tadpole is momentum independent the propagator can be parametrised as 
Gip{p) = if {p 2 — rnfy where is determined by the gap-equation 

2 2 , A3 2 , ^lm la \ 

m^ = m 2 + —v + —T^. (6) 

Finally, the closed system determining the propagator G a and the 3-point function Y^ a 
is derived by taking the corresponding second and third functional derivatives of the 
effective action: 



ST 

8a\ Vu =o 



Nv 



2 . ^2 2 A 3 

m 3 + —v +y T ^ 



2 6 

Aa 



Fw*(p,q, -p-q) 



-iy/Nv-^- / G^(p - k)G^(k)F^ a (p - k,k, -p), 

O Jk 

\ 3 V 



3VN 
.Ai 



-i— - / G^(p + q - k)G^(k)T^ a (p + q - k,k, -p - q). (7) 

Jk 



We note that T^ a was calculated from the master equation of ST/St/j. In the last equation 
the expression of the right hand side shows that T^ a (p, q, —p — q) on the left depends 
only on the momentum of a. From this it promptly follows for the 3-point function that 

IWrfo, g> -P-v) = I _ Mj^ p + ' I^(p) = -iJ k G^(p-k)G^(k), (8) 

which leads (by exploiting Eq.(4)) to 

■ri-it \ 2 A 2 2 A3 2 ^v(p) /n\ 
3 18 1 - ^(p) 



The renormalisation algorithm of the equations (4), (5) and (7) is rather non-trivial. The 
result can be stated, however, very simply: i) one should replace in these equations the bare 
couplings mf , \ by the corresponding renormalised ones; ii) the quadratically divergent 
tadpole integral and the logarithmically divergent 1^ should be replaced by the finite 
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parts T^^, Ii(i,f arising after subtracting the divergent pieces. It is worth noting, that by 
the renormalisation of \ the 3-point function T^ a is actually a finite asymptotically 
decreasing function of the a-momentum. As a consequence no divergence proportional to 
p 2 appears in G" 1 , therefore there is no need for infinite field renormalisation. 

The derivation of the non-trivial relation of the renormalised couplings to the unrenor- 
malised ones will be outlined in the next section. We apply cut-off regularisation to the 
divergent integrals T^(M) and I^(p, M), which are evaluated after Euclidean continuation 
of the variables: 



r i A 2 M 2 eA 2 

T * (M) s L w^w = is? - IS? ln MS + T * AM) ' 

f % 1 eA 2 

'*<"■ M ) - 1 ( B-KPnP-kf-M'') = -T^ ,n MS + '* Ap - m) < (10) 

where the scale M is introduced to eliminate the M-dependence of the cut-off dependent 
terms. The renormalisation will be achieved by compensating the cut-off dependence 
explicitly written out above by appropriately chosen coupling counterterms. 



3 Renormalisation 



The ultraviolet renormalisation of the equations (4), (5) and (7) requires a fresh analysis 
of the counter-term resummation, since the unrenormalised equations do not suggest any 
transparent formula similar to the simpler case of the O(N) model (see e.g. [13]). 

Below we shall apply the iterative renormalisation of [11] to our case. The renormali- 
sation of the equation of state and pion propagator (gap-equation) is the usual super- 
daisy renormalisation described in details in [11]. The renormalisation of the system (7) 
is a little bit more involved. First one rewrites the previous equations by making ex- 
plicit the counterterms. Since the derivation of the DS equations starts from the classi- 
cal action, there are explicit Lagrangian parameters present which we have to split as 
m 2 = m 2 R + 5m 2 , Aj = Aj^ + 5A«. But the effect of the counterterm does not reduce 
merely to this. Their consistent application implies that the 3-point function entering the 
DS equation of the sigma propagator is the renormalised one. This statement can be easily 
checked when using in the equation of the sigma propagator the iterative solution of the 
equation of the renormalised 3-point function: 



IW(p,g, -P-q) 



A3,i? + #A3 



— I 



3VN 
Ai,j? + 5Xi 
6 



11) 



/ Gip(p + q - k)G^(k)T ffa (p + q - k,k, -p - q). 

Jk 
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At a given order of the iteration we recover the chain of bubble diagrams including the 
counterterm diagrams of the usual perturbation theory at the corresponding order of the 
loop expansion. 

Next we subtract Eq. (4) from the first equation of (7) to obtain for the renormalised 
iG^ 1 the expression 

*G-\ P ) = P 2 - X2 ' R ± 5X \ 2 -iyfNv Xs ' R ± 6X3 I G4p-k)G4k)T^ a (p-k,k,-p). (12) 

O J k 

The iterative renormalisation consists of assuming infinite series expansion for the coun- 
terterms and determining them by requiring the cancellation of the divergencies generated 
at each iteration step of the solution of the coupled equations of G a and Y^ a . 

In the first iteration of the vertex function one substitutes into the propagator equation 
the tree-level value of the 3-point vertex which yields the requirement: 

+ %V (-< / G^P - k)G^k)) = 0. (13) 

O 10 V Jk J div 

The same sort of first iteration is done also in the integral equation of T^ a and the two 
together lead to 

= ^r-I^div, = 1,R 3 ' R I^,div (14) 

With this step also the finite parts are determined, which are substituted into the second 
iteration. This step produces on the right hand side of the expression iG^ 1 the following 
divergent piece in the full expression, which has to vanish: 



. r,*3,Ru*3 t , *3,R T i s I n (l) . *l,R*3,R T \ . ^1,R^3,R T 2 
+ Z h,div + —^-h,FK.P) \0*3 H g H,div I H J ^,div 



v 2 = 0. 



dX^ , ^X^rSX^ t t A 3;jR T f^fx^ii) , Ai )R A 3jR T \ t X 1>R \l 

(15) 

It is important that the subdivergence proportional to the finite (possibly temperature 
or density dependent) I^^p automatically vanishes in view of the previous iteration step. 
This feature of subdivergence cancellation is generally valid for each iteration step. The 
remaining divergent terms determine the second term in the infinite series of <5A 2 : 

sx?) = ^£, idn , (16) 

The second iteration of the equation of determines 5X^ and 8X^ = —X\ R I^^iv/^- 
With simple induction the following general recursions are found for the subsequent iter- 
ative terms of the counterterm series: 

JT\(") _ ^ n ~ l )^ R T A\( n > - X\(n-1)X 3 ,R ..( n ) _ ,.(n-l)A 3 ,il 

0*1 — —0A 1 — — l^div, — ~ 0A 3 —^--liprfiv, 0A 3 ~ —0A\ —g— 

(17) 
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where n > 1 and 5X^ := Ai ; #, := A 3v r. 

These equations imply the following non-perturbative renormalisation formulae: 



1 1 1 

Ai ~ \!, R + 6 



A3 
Ai 



Ai. 



A, 



A '-A 

a7 



2,R 



X 2 



Al 



(18) 



In addition to these relations between bare and renormalised couplings we also obtain the 
renormalised sigma self-energy which is nothing but the expansion of the last term in (9) 
in powers of Ai up to the given order of the iteration (of course, the couplings and the 
bubble integral have to be replaced by their renormalised expressions). 

It is easy to check that provided we would have known a priori these relations, they indeed 
ensure the simple form of the renormalised expression of (9) announced at the end of the 
previous section. The mass of the Higgs particle is determined by the pole of equation (9) 
(*G~V = M 2 )=0) 



Mi = — 

a 3 



A 



2,R 



X 2 s 
A 3,R 



Al. 



+ 



X 2 

A 3,R 



R, 



Kr K'r ~ hAP 2 = M2)/G\ 



(19) 



The renormalised quartic couplings are used also in the iterative renormalisation of the 
equations for the order parameter v and of in which the tadpole is calculated with a 
self-consistent propagator. The mass counterterms b~m\ an d $ m 3 are determined as: 



6X1 



(n-l) 



A 2 



5X 



(n-l) 



167T 2 

A 2 



+ m 2,R^ 
A 



SX^ 



(n-l) 



A 2 



3,R 



6 16tt 2 



2 SX\ 

+ ™2,R 



(n) 



6 167T 2 



+ m 2,R- 



SX^ 



A. 



3,R 



8Xt^> A 2 

6 167T 2 



Xl,R 

(")■ 



2 ^1 ^3,fl 
+ m 2,R x2 

A 1,R 



(20) 



implying 



m| A 2 
A7 + 96^ 



m 



2,R 



2 2^3 2 2 ^3,R / 01 \ 

m 3~ m 2Y- = m 3,R ~ m 2,R-, • ( 21 ) 

M,R A l M,R 

The second equalities in both equations of (20) follow when one exploits (17). Making 
use of these renormalised relations in (4) and (5) one easily derives for the renormalised 
expression of the vacuum expectation value the following formula: 



A 



— m 



3,R 



+ [ m 2,R ~ m 



A 



3,R 



Ai 



,R 



X2,R — 



X 2 x 
A 3,R 



Al 



■R, 



(22) 



4 Discussion 



The expressions for the Higgs mass (19) and the Higgs condensate (22) show that the 
influence of the phantom sector might be quite strong on both. In order to make the 
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discussion more transparent below we shall choose 



A 



2,R 



^3,R = K ^1,R = A + A'. 



(23) 



In this case the gap equation (6) which determines reads after the elimination of v 2 

as 



m 



m 2,R m 3,R 



A'. 

6~ 



-T, 



T, 



rn 



4>,F 



TP 



16tt 2 



In 



which has the approximate solution (for small A') 



em 



1> 



Mi 



(24) 



m 2,R m 3,R 



A' j n < m l,R- m l,R) ' 



1 - 



96vr 2 



(25) 



For sufficiently large renormalisation scale M one has < m 2 < m\ R — m\ R . The Higgs 
condensate has the simple equation 



^2 2 ^ / 2 

-v = -m 3R + —(m 2R 



m 



3,R 



ml). 



(26) 



Depending on the ratio of the two quartic couplings and the differences in the renor- 
malised mass parameters we see that the phantom fluctuations might produce a stable Z 2 
symmetry breaking condensate. This mechanism is similar to the dynamical symmetry 
breaking of Coleman and Weinberg [14]. The role of the electromagnetic field in their 
example is played here by ipi, A' and the dynamically determined m 2 . Also the Higgs 
mass 



Ml 



A 



3 A + A' 



A' + 



A 



A+A 



-I^ F (P 2 = Ml) 



(27) 



depends in a substantial way on both couplings which would give more flexibility in Higgs 
phenomenology. 



In conclusion, we have presented an O(N) symmetric model of the phantom world in- 
teracting with the SM Higgs field, which was solved in the N = oo limit. The solution 
obtained with iterative renormalisation represents a non-trivial generalisation of the well- 
known example of the pure O(N) model. The method can be applied to a large variety 
of models containing several multiplets with large number of components. 

The existence of symmetry breaking solutions stabilised by the phantom sector fluctu- 
ations was pointed out. In this respect we emphasise that Eqs. (24), (26) and (27) are 
valid also at finite temperature, if I^p and T^p are evaluated accordingly. They can serve 
for investigating the nature of finite temperature symmetry restoration in the generalised 
Higgs sector. 

Non-perturbative (lattice) studies at finite iV could clarify if our findings are generically 
valid also for finite number of phantom field components. 
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